Partonic Picture of Nuclear Shadowing at Small x 
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We investigate the nuclear shadowing mechanism in the context of perturbative QCD and the 
Glauber multiple scattering model. Using recent HERA data on nucleon structure function at small 
X, we put stringent constrains on the nucleon gluon density in the double-logarithm approximation. 
We suggest that the scaling violation of the nucleon structure function in the region of small x and 
semihard scale can be reliably described by perturbative QCD which is a central key to the 
understanding of the scale dependence of the nuclear shadowing effect. Our results indicate that 
while the shadowing of the quark density arises from an interplay between the "soft" and semihard 
QCD processes, the gluon shadowing is largely driven by a perturbative shadowing mechanism. We 
demonstrate that the gluon shadowing is a robust phenomenon at large and can be unambiguously 
predicted by perturbative QCD. 
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I. INTRODUCTION 

The theoretical studies of semihard processes play a 
increasingly important role in ultrarelativistic heavy-ion 
collisions at ^/s > 200 GeV in describing the global colli- 
sion features such as the particle multiplicities and trans- 
verse energy distributions. These semihard processes can 
be reliably calculated in the framework of perturbative 
QCD and they are crucial for determining the initial con- 
dition for the possible formation of a quark-gluon plasma. 
Assuming the validity of the factorization theorem in per- 
turbation theory, it is essential to know the parton distri- 
butions in nuclei in order to compute these processes. For 
example, to calculate the minijet production for 2 
GeV or the charm quark production in the central rapid- 
ity region (y = 0), one needs to know the nuclear gluon 
structure function, xgA(x,Q^), &t Q ^ pT or wic and 
X = x-B] = 2pTl\fs or 2mcl\fs- For ^/s > 200 GeV, 
the corresponding value of a; is a; < 10~^. Similarly, the 
knowledge of nuclear quark density xfA{x, Q^) is needed 
for the Drell-Yan Icpton-pair production process. 

The attenuation of quark density has been firmly es- 
tablished experimentally in deeply-inelastic lepton scat- 
terings (DIS) from nuclei at CERN ^ and Fermilab 
1^ in the region of small values of the Bjorken variable 
X = Q^/2mv , where is the four-momentum trans- 
fer squared, v is the energy loss of the lepton in the 
lab frame and m is the nucleon mass. The data, taken 
over a wide kinematic range 10~^ < a; < 0.1 and 0.05 
GeV^ < < 100 GeV^, show a systematic reduction 
of nuclear structure function F^{x,Q'^)/A with respect 
to the free nucleon structure function {x, Q^)- There 
are also some indications of nuclear gluon shadowing from 
the analysis of J/^* suppression in hadron- nucleus experi- 
ments j^] . Unfortunately, the extraction of nuclear gluon 
density is not unambiguous since it involves the evalu- 
ation of initial parton energy loss and final state inter- 
actions Q. These measurements present a tremendous 
challenge to the theoretical models of nuclear shadowing 



(for recent reviews, see 0). 

The partonic description of nuclear shadowing phe- 
nomenon has been extensively investigated in literatures 
in the context of a recombination model and a 

Glauber multiple scattering theory |§-|l^. These ap- 
proaches predict a modified QCD evolution equation for 
the nuclear structure function in contrast to a Gribov- 
Lipatov-Altarelli-Parisi (CLAP) evolution equation |lj] 
for nucleon. Qualitatively, the scale (Q^) evolution of 
the nuclear structure function in the modified equation 
is slower than the nucleon case, leading to a perturba- 
tive mechanism for the depletion of the parton density 
at large momentum transfer. The quantitative predic- 
tion of nuclear shadowing at large in general de- 
pends on the initial value of shadowing at small mo- 
mentum transfer where perturbative calculation breaks 
down and some non-perturbative model or experimental 
input is needed. The so-called vector-meson-dominance 
(VMD) model ||T^ has been quite successful in predicting 
the quark shadowing at low where the experimental 
data are rich ||l|,|^. Presently, there is no viable non- 
perturbative model for the gluon case and there is very 
little experimental information on the initial gluon shad- 
owing value at some low Qq- This, in practice, would 
make the prediction on gluon shadowing very uncertain. 

In this paper, however, we suggest that there exists a 
typical scale ~ 3 GeV^ in the perturbative evolution 
of the nuclear gluon density beyond which the nuclear 
gluon shadowing can be unambiguously predicted in the 
context of perturbative QCD. The so-called "semihard" 
scale Qgjj is determined by the strength of the scaling vi- 
olation of the nucleon gluon density in the small- a; region 
at which the anomalous dimension, 7 = dlnxgj^ / XnQ^ , 
is of 0(1). We shall demonstrate that as approaches 
Qgjj, the gluon shadowing ratio rapidly approaches a 
unique perturbative value determined by the Glauber 
multiple scattering theory. For > QgH' the nuclear 
gluon density, in contrast to the quark case, is almost 
independent of the initial distribution at Qq and its evo- 
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lution is mainly governed by a GLAP equation. The 
suggested phenomenon is a consequence of the singular 
growth of the gluon density in the small-x region and the 
perturbative shadowing caused by the coherent parton 
multiple scattering. 

In particular, we shall present a quantitative study 
of the scale dependence of the nuclear shadowing phe- 
nomenon for both the quark and the gluon distributions 
based on the Glauber diffractive approximation. A simi- 
lar study has been performed by Eskola in the recombina- 
tion model . Our focus is the sensitivity of the nuclear 
shadowing at large to the initial conditions. To de- 
termine the gluon content of the nucleon, we shall use 
a leading double logarithm approximation (DLA) where 
the scaling violation of the quark structure function 
is entirely caused by a g — > splitting. The agreement 
between our results based on the DLA and the measured 
values of indicates that the perturbative evolution 
sets in as early as > 0.8 GeV^ and the scaling vio- 
lation of gluon structure function in the semihard region 
can be reliably described by the DLA. This provides a 
basis for studying both the nuclear quark and gluon den- 
sities in the Glauber approach where the separation of 
qq or gg pair in the DLA is considered frozen during its 
passage in the nucleus. The reason for not using the 
available parton distributions for a nucleon such as those 
of MRS, CTEQ or GRV parameterization ^ is that the 
nuclear structure function obtained in our Glauber ap- 
proach does not contain a full GLAP kernel, as it will 
become clear below. Since the shadowing ratio depends 
on both the nuclear and the nucleon structure functions, 
they should be consistently calculated in the same ap- 
proximation. We expect that the calculated ratio, how- 
ever, has a validity beyond the DLA. 

We shall demonstrate that the nuclear shadowing for 
the quark distribution is a leading twist effect due to the 
interplay between the "soft" and the "semihard" QCD 
processes, while for the gluon distribution, the shadowing 
is largely driven by the perturbative shadowing mecha- 
nism at small x. Beyond the semihard scale, both quark 
and gluon shadowing ratios are shown to evolve slowly. 
We discuss the x-dependence of the shadowing ratio and 
the non-saturating feature of the perturbative shadowing. 
As a by-product of the Glauber approach, we also study 
the nuclear geometry dependence of the shadowing ratio, 
which cannot be addressed in the recombination model 
such as the GLR equation. Since the shadowing effect is 
intimately related to the effective number of scatterings, 
it should vary strongly as one changes the impact param- 
eter. This leads to the breakdown of the factorization of 
the nuclear parton density into the product of an impact 
parameter averaged parton distribution and the nuclear 
thickness function. 

The paper is organized as follows. In section ^ we dis- 
cuss the physical picture of the nuclear shadowing phe- 
nomenon. It can be regarded as a heuristic introduction 



to the central idea of our approach. In section III, we cal- 



culate the nucleon structure function F^{x, Q^) for large 
range of x and using a gluon density obtained in the 
DLA. We determine the initial condition for and for 
xgN by fitting the Q^-depcndcnce of F^ at small value 
of X. The gluon density in a nucleon is crucial to the 
nuclear shadowing effect. In section IV we give a brief 



introduction on the Glauber multiple scattering model 
and discuss the dependence of the nucleon density on 
the coherence length of the virtual photon. The nuclear 
shadowing effect for F^ is calculated in section ^ and 
the initial shadowing is parameterized in order to fit the 
available experimental data at the measured values. 
In section VI, the gluon shadowing is computed using 
the Glauber formula derived by Muell er and Ayala 
F, Gay Ducati and Levin p^. Section VII is devoted to 



the discussion of the impact parameter dependence of the 
nuclear shadowing. We conclude our findings in section 
VIII. 



II. PHYSICAL PICTURE OF NUCLEAR 
SHADOWING 

At low {< 1 GeV^) in the deeply-inelastic scat- 
tering, the interaction of the virtual photon with the 
nucleons in the rest frame of the target is most nat- 
urally described by a vector-meson-dominance (VMD) 
model psf , in which the virtual photon interacts with 
the nucleons via its hadronic fluctuations, namely the 
p, Lu and (j) mesons. The shadowing of the (sea) quark 
density is caused by the coherent multiple scatterings of 
the vector meson off the nucleons with a large hadronic 
cross section avN, which can be most conveniently stud- 
ied in the Glauber diffractive approximation model [ p^ . 
The x-dependence of the shadowing ratio, R{x,Q^) = 
F^{x,Q'^)/AF.^ {x,Q^), is explained by the coherent 
length of the virtual photon Ic — l/2mx. For x < 10^^, 
Ic > Ra = 4 ^ 5 fm for large nuclei, the vector me- 
son interacts coherently with all nucleons inside the nu- 
cleus resulting in a uniform (in x) reduction of the total 
cross section, while for x > 10"^, Ic < Ra the vector 
meson can only interact coherently with a fraction of 
all the nucleons and shadowing is reduced. The VMD 
model in general predicts much less scaling violation of 
F^{x,Q'^) than it has been observed in the small x re- 
gion and a rapid disappearance of the nuclear shadow- 
ing for ^ niy. Although a generalized VMD model 
which includes higher resonances may improve the situ- 
ation somewhat, it is generally believed that the nuclear 
shadowing in the VMD model is a higher twist effect, 
which is in contraction with the experimental evidence 
that the R{x, Q^) has a rather weak dependence 

At high (> 1 2 GeV^), the virtual photon can 
penetrate the nucleon and probe the partonic degrees of 
freedom where a partonic interpretation based on pertur- 
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bative QCD is most relevant in the infinite momentum 
frame. It is well known that the Q^-dependence, the 
so-called scaling violation of the nucleon structure func- 
tion {x, Q^) can be well accounted for by the Gribov- 
Lipatov-Altarelli-Parisi (GLAP) evolution equation flj] 
given some non-perturbative initial condition F^{x,Q^) 
[ pTf . In the small x region, the GLAP equation pre- 
dicts a strong scaling violation of {x, Q^) due to the 
strong growth of the gluon density in the nucleon. In 
fact, a dominant perturbative evolution implies that the 
growth of F^{x,Q'^) in InQ^ at fixed, small x is closely 
related to the growth in Inl/x at fixed Q^, i.e. the so- 
called double scaling violation observed in the DIS ex- 
periments for nucleon, most transparent in the Leading- 
Double-Logarithm Approximation (DLA) of the gluon 
structure function. At extremely small x, the rcsum- 
mation of large \nl/x is necessary leading to the BFKL 
evolution equation ||l9|] which may also be derived in a 
Weiszacker- Williams classical field approximation ||2^ . It 
is not clear when the BFKL hard pomeron will become 
relevant since the GLAP evolution seems able to describe 
most of kinematical region that current experiments have 
probed so far. As we shall show, the main feature of the 
Q^-dependence of the nuclear shadowing is intimately re- 
lated to this strong scaling violation in the small-x region. 

If the nuclear structure function follows the same per- 
turbative GLAP evolution, given an initially shadowed 
F^{x, Qq), the shadowing ratio R{x, Q^) as increases 
will quickly approach 1 because of the strong scaling 
violation F^{x,Q'^) > F^{x,Ql). This would imply 
that the shadowing is purely non-perturbative and thus 
a higher twist effect. The key observation is that in the 
infinite momentum frame, because of the high nucleon 
and parton densities, the quarks and gluons that belong 
to different nucleons in the nucleus will recombine and 
annihilate, leading to the so-called recombination effect 
first suggested by Gribov, Levin and Ryskin [|| and later 
proven by Mueller and Qiu 0. In the target rest frame, 
the virtual photon interacts with nucleons via its quark- 
antiquark pair (qq) color-singlet fiuctuation [|lo|. If the 
coherence length of the virtual photon is larger than the 
radius of the nucleus, Ic > Ra, the qq configuration in- 
teracts coherently with all nucleon with a cross section 
given by the color transparency mechanism for a point 
like color-singlet configuration ||2l[]. That is, the cross sec- 
tion is proportional to the transverse separation squared. 



<JqqA - / d'^b2{l - . 
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(2.2) 



where b is the impact parameter and TA{b) is the nuclear 
thickness function. The nuclear cross section is therefore 
reduced when compared to the simple addition of free 
nucleon cross sections. 



However, it is not clear that (2.1) and (2.2) will provide 
a satisfactory explanation of the non-vanishing shadow- 
ing effect at large Q^. Although for a given Q^, all qq 
configurations with the separation size > are 
possible, the virtual photon wavefunction |V'7*P ^ V^t 
favors the smallest value of ~ 1 /Q^ where the cross 
section is small and no significant shadowing would be 
expected. The resolution to this problem comes from the 
fact that the nucleon structure function exhibits a strong 
scaling violation in the small x region at the moderate 
(semihard) momentum transfer 0.8 GeV^ < < Qgn 
where Qg^ ~ 3 GeV^. As we shall show, the strong 
scaling violation in this region can be understood in the 
context of perturbative QCD. The gluon density in (2.1) 



grows rapidly with l/r| in the small x region, which can 
compensate the inherent "smallness" of the cross sections 
of "hard" QCD processes. As a consequence, the semi- 
hard processes can compete successfully with the "soft" 
processes as described in the VMD model and continue 
to provide the shadowing mechanism up to some moder- 
ately high value of ^ Qsa ■ terms of the anomalous 
dimension 7(<5^) = dlnxgN /dlnQ'^ , the cross section in 
may be effective expressed as 



aqqN OC 



(Q 



2\j 



(2.3) 



The large anomalous dimension, 7 > 1, even at moder- 
ately large delays the rapid falling of the cross section 
due to the "smallness" of the configuration size, as first 
argued by Gribov, Levin and Ryskin in We refer to 
this as the perturbative shadowing mechanism as it can 
be computed in perturbative QCD. At even higher val- 
ues of Q^, > QgHi the scaling violation becomes weak 
(logarithmic), and the perturbative shadowing mecha- 
nism ceases to be effective. In this case, both the nucleon 
and the nuclear structure function evolve slowly at the 
same rate. The shadowing ratio developed by both per- 
turbative and non-perturbative mechanisms in the rela- 
tively low-Q^ region persists in the high-Q^ region but 



of the q and q. In the double logarithm approxima- tends to 1 rather slowly as Q 00. 



tion (DLA), the cross section can be expressed through 
the gluon structure function for a nucleon 



cTqqN OC Usx' guLAix , l/rj) , 



(2.1) 



where x' = K'P /2nnv and M is the invariant mass of qq 
pair. The total cross section GqqA can be calculated in 
the Glauber eikonal approximation. 



In contrast to the VMD model, the x-dependent cross 
section ciqqiq predicts an additional x-dependence of the 
shadowing ratio apart from the coherence length depen- 
dence in the perturbative region. As a result, one does 
not in general expect the saturation of the shadowing ra- 
tio as X decreases below the coherence limit x < 10^^. 
On the other hand, as decreases, the gluon distri- 
bution changes from a singular small-x behavior (hard 
pomeron regime) to the less singular one (soft pomeron 
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regime). The shadowing ratio becomes flatter as de- 
creases. The observed saturation of the shadowing ratio 
in the small- a: region Q is thus a non-perturbative, low- 
feature. 

Similarly, the nuclear gluon density may be studied by 
considering a deeply-inelastic scattering of a virtual col- 
orless probe (e.g. the virtual "Higgs boson") with the 
nucleus where the interaction proceeds via a gluon pair 
[gg) component of the virtual probe. This has been in- 
troduced by Mueller |jl^ and recently studied by Ayala 
F, Gay Ducati and Levin (AGL) In contrast to 

the quark distribution, the gluon structure functions for 
nucleon and nucleus are poorly known, especially in the 
small- cc region. They are both of extreme importance 
for the semihard QCD processes at high energies. Most 
of our knowledge about the gluon distribution comes 
from the theoretical expectations based on perturbative 
QCD constrained by the scaling violation of F2 and the 
prompt photon experiments . The nuclear gluon den- 
sity may also be studied in the Glauber multiple scatter- 
ing model where Oqqjss in (2^) is to be replaced by Oggjq. 
Simple color group representation argument relates these 
two cross sections by Uggfq = 'duqqt^j^. However, what 
makes the gluon case very different from the quark case 
is the unknown initial (non-perturbative) gluon density 
in the nucleus. There is very little experimental informa- 
tion on the gluon shadowing at low and small x. The 
VMD model is not applicable here since the low-lying me- 
son spectroscopy for gg-\^k& states is poorly known. This 
would make this approach very uncertain in predicting 
the amount of gluon shadowing in the nucleus. 

However, as we shall show in this paper, the larger 
cross section cTg^jv implies an even stronger scaling vio- 
lation of gluon density than the quark case, making the 
perturbative shadowing a dominant mechanism. As a 
consequence, the shadowing ratio is largely determined 
by the perturbative calculation at — Qgjj- That 
is, the large anomalous dimension of gluon distribution 
function for a nucleon drives the shadowing ratio rapidly 
to approach the perturbative value at moderately large 
— Qgfj the influence of an initial shadowing at low 
becomes relatively unimportant. For > QIh' 
scaling violation becomes weaker, the shadowing ratio 
determined by the perturbative shadowing mechanism 
at the semihard scale persists and evolves slowly. The x- 
dependence of the shadowing ratio can be also predicted 
to be non-saturating, a generic feature of the perturba- 
tive shadowing. We therefore suggest that the nuclear 
gluon shadowing at > Qgjj ~ 3 GeV^ and at small x 
can be predicted in the context of perturbative QCD. The 
prediction is based on the scaling violation of the gluon 
structure function and the existence of the perturbative 
shadowing in the Glauber multiple scattering model. 



III. LAB FRAME DESCRIPTION OF NUCLEON 
STRUCTURE FUNCTION 

Since the nuclear shadowing ratio depends sensitively 
on the nucleon structure function, we shall develop a for- 
malism to calculate both the nuclear and the nucleon 
structure functions simultaneously in order to reduce the 
uncertainty due to the nucleon structure function. This 
prevents us from using the well-known parametrizations 
such as MRS, CTEQ and GRV which are based on 
the full next-to-leading-order CLAP equation because 
the nuclear structure function is not derived in the same 
formalism. We stress that the shadowing is only a few 
tens percentage effect, thus the uncertainty in the nu- 
cleon structure function at the same level can cause a 
large effect on the shadowing ratio. 



A. Quark Distribution 

Consider a description of deeply-inelastic lepton scat- 
tering off a nucleon in the laboratory frame. The inclu- 
sive cross section can be expressed via the virtual photon- 
nucleon total cross section. The cross section for the 
absorption of a virtual photon in the small x region is 
dominated by the scattering off the sea-quark qq pair. 
The generic perturbative QCD diagrams giving rise to 
the qq fluctuation are shown in Fig. |l|. The CM energy 
of the incoming virtual photon at small x is 



s = (q+p)^ :^2p- q 



(3.1) 



where q and p are the four-momenta of the photon and 
the target nucleon. Thus a region of small x corresponds 
to a high energy scattering process at fixed Q^. 





FIG. 1. The perturbative diagrams giving rise to the qq 
pair fluctuation. 

The imaginary part of the amplitudes in Fig. |l| is re- 
lated to the photoabsorption cross section which has been 
calculated by Nikolaev and Zakharov assuming that 
the size of the qq pair is frozen in the scattering pro- 
cess and that the one-gluon exchange process dominates. 
This is nothing but the so-called double logarithm ap- 
proximation (DLA) The transverse cross section 
can be cast into an impact parameter representation 



a{j*N) = / dz 
Jo 



(fr\'tl;{z,r)\^aggN{r) 



(3.2) 
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where z is the Sudakov variable, defined to be the frac- 
tion of the qq pair momentum carried by one of the 
(anti)quarks, r is the variable conjugate to the trans- 
verse momentum of (anti)quark of the pair, representing 
the transverse size of the pair, ■(/'(•^^j r) can be interpreted 
as the wave function of the qq fluctuations of virtual pho- 
ton. Thus, I'l/'C^jr)!^ serves as the probability of finding 
a qq pair with a separation r and a fractional momentum 
z. In an explicit form 



efP{z)a^Ki{ary 



(3.3) 



where = Q^z{l — z) and Ki is the modified Bessel 
function. The factor P{z) = + (1 — z)^ is the universal 
splitting function of gauge boson into the fermion pair. 
Since Ki is a rapidly falling function, the integral in ( |3.2| ) 
has important contributions only in the region ar ^ I. 
In this case, one has 



Nf 

Oaem 2 J_ 

(27r)2 2^^* r2 



(3.4) 



where the expansion Ki^ar)"^ has been used. 

The photon wave function favors a qq pair with small 
transverse separation. 

The cross section for the high energy interaction of a 
small-size qq configuration with the nuclcon can be un- 
ambiguously calculated in QCD in the region of small x 
by applying the QCD factorization theorem. In the DLA, 
the result is [E3| 



where xg^^^ is the gluon distribution of the nucleon cal- 
culated in the DLA. The fact that the cross section is 
proportional to the size of the point- like configuration is 
the consequence of the color transparency in QCD. Al- 
though the factor in (J^^ decreases fast as one goes 
to the short distance, the singular behavior of the wave 
function in (3^) and the strong scaling violation of xg^^^ 
in the small- a; region as r decreases can compensate the 
smallness of cross section due to the color transparency. 
As a result, the soft and hard physics compete in a wide 
range of Q^. 

The x' and Q'^ in (3.5) are related to z, x and as 
follows. The invariant mass squared of the qq pair is 



z(l-z) ' 



(3.6) 



where fcj ~ is the transverse momentum squared 

of the (anti)quark. The virtuality of the gluon is then 
Q'^ = {kt — {—kt)Y = while x' is given by 



M 



-f2 



where u is the energy loss of the lepton. The condition 
ar < 1 implies x' > x. Since > 4/(5^, the integration 
region in z is limited by a^r^ = Q'^r^z{\ — z) < 1. Be- 



cause the integrand in (3.2) is symmetric under z <J=> 1 — z, 
one has dz ^ 2 dz where e is the solution to ar = 1 
or x' — X. 

At small X, the free nucleon structure function can 
be written in terms of the virtual photon-nucleon cross 
section 



Fi'ix, 



efx/Nix^Q'^ 



(3.8) 



where xf^ix^Q^) is the (sea) quark distribution func- 
tion. Substituting (3.2) into the above equation and 
changing the integration variable z into x' in (3.2) dz = 
—xdx'/{x''^Q^r^)^ one obtains 

fN{x,Q^) 



^ , 3 r dx' 



4/Q2 



-(^qqN{x',r ) 



(3.9) 



where aqqN is given in ([3.51) . We have split the integration 
on into the region 4/Q^ < < 4/(5§ where the per- 
turbative calculation is valid and the region > A/Qq 
where the perturbation theory breaks down and the cross 
section is not calculable. The non-perturbative part is 
taken as an initial condition /Ar(a;,Qo) which should be 
determined by the experimental measurement. However, 
it is worthwhile to emphasize that because of the expo- 
nential cutoff of the large by Ki{ar)^ in (3.2), the inte- 
gration in is infrared safe. Even if one substitutes the 
small-r cross section into ( |3.9| ), the integral is only log- 
arithmically divergent ~ It^Qq, in contrast to the usual 
factorization scheme of the hard process. This allows a 
smooth extrapolation of the "semihard" physics into the 
"soft" physics. Thus one can choose a rather small cutoff 
Qq in order to integrate more perturbative contributions. 
Because of the infrared safety, the uncertainty due to the 
choice of Qq is minimal. 

Substituting (^) into ( ^!9| ) and replacing 4/r^ by Q'^, 
one recovers the familiar integrated form of the CLAP 
evolution equation at small x 



1 



fN{x,Q')^fNix,Q'Q) + — / — 



dx' 



dQ'- 
Q'2 



■«.(Q")ff^^^(^',Q"). (3.10) 



Note that ( 3.1Cl| ) is not an integral equation since it does 
not involve the sea-quark contribution on the right hand 
side. It is derived in the DLA where the small- cc behavior 
is driven by the gluon contribution. To improve the larger 



2mv z(l — z)2m^ a?r'^ 



(3.7) 

X behavior, we insert a splitting function Pgq(z) — 
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(1 — z)^ where z = x/x' in ( |3.10 ) and add the valence 
quark contributions to F2 



l/x. To improve the DLA at a large (> 20 GeV^), 
the resummation of single logarithm, as In Q"^ IQ\, is nec- 
essary. This results in an additional dependence of 



F^{x, Q^) = ^ ef[xfN{x, Q^) + xvi{x, Q^] , (3.11) the gluon distribution, a multiplicative factor (t/to) 



where xv^s are the valence quark distributions, which 
are given by the GRV |Q parameterization (there is a 
very little difference among different parameterizations 
for valence quarks). 



( p^ where 6 = 61/45 We find that this mod- 

ification has very little numerical effect on the nuclear 
shadowing ratio in the region of interest. 

In a semiclassical form xqn — a;~'^(Q^)^. One defines 
the exponents 7 = dlnxQN /dt and uj — dxgN / dh\{l/ x). 
In the DLA, 



B. Scaling Violation of F^{x,Q^) and Gluon Density 

The evolution or the scaling violation of the nuclcon 
gluon density in the small-x region is quite different from 
that in the large- a; region in the context of GLAP equa- 
tion. In the region of small x, one has to sum up large 
tts In 1/a; contributions in addition to as In terms rep- 
resented by the ladder diagrams with strong scale or- 
dering in the leading logarithm approximation (LLA). 
To avoid summing up two large logarithms, one consid- 
ers the DLA, i.e., the most leading diagrams involving 
(a^ In In l/x)" at each order n. Such a double loga- 
rithm appears only when vector particles (gluons) propa- 
gate in the i-channel of the ladder diagram as argued by 
Gribov, Levin and Ryskin in Thus, the scale depen- 
dence of the nucleon quark distribution is directly related 
the small- a; behavior of the gluon distribution. 

In general, the DLA is valid in the kinematical region 
where ofg Inl/xlnQ^/Qg ~ 1 while In l/x <C 1 and 
as IhQ^/Qq <C 1 as well as as <C 1. The precise validity 
range of DLA or of GLAP equation is currently still under 
debate. The general belief is that it describes a region up 
to a moderately large (^ 20 GeV^) and a moderately 
small X, 10~^ < X < 10~^. Formally, in the small-x 
region, one should sum the large as In l/x terms using 
a BFKL resummation technique |19|. However, it has 
been shown recently by Ball and Forte ||2^ that a gluon 
density obtained in the DLA together with the resummed 
leading logarithm in gives a very good description of 
the double-asymptotic scaling observed at HERA [ p4[ . 
The unitarized gluon density in the DLA has a simple 
approximate form P,p2[, 



DLA 

xgN 



xgN{x,Ql)Io{y) , 



(3.12) 



where /q is the Bessel function and 



with t = \nQ^ I and c = ^36/25, provided that input 
xg^ix, Qq) is not singular in x. The asymptotic behavior 
of Jo(y) is that Io{y) — > 1 as y — > and Io{y) j ^^-Ky 
as y 00. That is, in this DLA, xgN increases faster 
than any power of In l/x, but slower than a power of 
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In l/x 
tint /to 



OJ 



(3.14) 



The anomalous dimension can be large in the small-x 
and -t region. As t ^ tg, a; — > and 7 is finite. It is 
important to emphasize that the DLA form (3.12) is only 
vahd when » Qq- 

Our goal is to pin down the gluon density in the 



nucleon. In our approach, the slope of K 



N 



dF^ /dlnQ^ , is virtually a direct indicator of the magni- 
tude of the gluon distribution. We require that the gluon 
distribution in the DLA as given in (3.12) describes rea- 
sonably well the scaling violation of for ranging 
from 0.85 GeV^ to 100 GeV^ where the experimental 
data exist in the small-x region. As it is turned out, this 
almost uniquely determines the gluon distribution. 
There are two unknown quantities in ( 3.12| ): Qq and 



the initial condition xgpfixjQo)- Qq determines where 
the perturbative calculation should begin to be valid and 
where the non-perturbative contribution can be sepa- 
rated into an initial condition x5Ar(x, Qg)- We would 
like to have a reliable perturbative description of the 
structure function in the range from 1 to a few GeV 
which is extremely important in understanding the evo- 
lution of nuclear shadowing. This leads to choosing 
Qo to be rather small so that the DLA sets in at a 
relatively low = 1 ~ 2 GeV^. Since xgN{x,Ql) 
has to be non-singular in the small x regime, we take 
x(7jv(x, Qq) = CgXg'^^^ {x, Qo) at an initial scale Qq since 
the initial gluon distribution in GRV is fairly flat while it 
has the proper behavior at x — > 1. The initial condition 
for is taken to be in the form motivated by a soft 
pomeron with an intercept a = 1 + 0.08 



Fi'{x,Ql) = cpx- 



^(l-x)" 



(3.15) 



where 77 = 10 accounts for the behavior as x ^ 1. We 
determine Qq, Cg and cp by requiring that the calculated 



(3.13) Q -dependence of at x = 8 



X 10~^ fits the recent 
measurement from HI Gollaboration |24|. We find the 
following unique solution 



ro = 0.4 ; 



Cg = 2.2 



CP = 0.2 



(3.16) 



The Q^-dependence of F^ at other values of small x are 
then predicted using ( 3.1C ) . In Fig. ||, we show our results 
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for 



compared to the experimental values The 
gluon density in ( 3.12| ) describes the scale dependence of 
F2 very well starting at a relatively low Q^, Q"^ ~ 0.85 
GeV2. 




10 

Q' (GeV') 

FIG. 2. The Q^-dependence of the nucleon structure func- 
tion calculated using the DLA and compared with experiment 
data [24]. 

The x-dependence of the nucleon structure function 
can also be predicted from ( 3.1Cl| ). Fig. ^ shows the cal- 



culated F.^{x,Q^) from (3.10) as function of x at var- 



ious as compared to HI measured values |24 . The 



agreement is evident. Starting with a fairly flat distribu- 
tion F.^{x,Qq), the nucleon structure function develops 
a more and more singular shape in x from the pcrturba- 
tive evolution governed by the gluon distribution. Some 
deviation is visible for x > 10~^, a region where we do 
not expect DLA to be valid. 



: Q'=0.85GeV' = 
- -- DLA gluon - 

E 


: 1.2 = 


1.5 


Jllllllll lllllllll lllllllll III 
I 3.5 : 


Jllllllll lllllllll lllllllll III 

; 5 : 

E 


- 1 iiiiiiii 1 iiiiiii 1 iiiiiii 1 li 
; 8.5 ; 


Jllllllll lllllllll lllllllll III 

: 15 : 

: 


Jllllllll lllllllll lllllllll III 
; 20 : 


- 1 iiiiiiii 1 iiiiiiii 1 iiiiiiii 1 li 

25 

■■ : 





XXX 

FIG. 3. The a;-dependence of the nucleon structure func- 
tion calculated using the DLA and compared with the exper- 
iment data. 



The results shown in Fig. ^and Fig. || demonstrate that 
the gluon density in the DLA given in (3.12) faithfully 
captures the essential feature of Q^- and a:-evolution of 
the nucleon structure function. The fact that we can 
describe even at very low Q"^, 0.85 < < 5 GeV^ 
indicates that the perturbative approach is reliable in the 
semihard region. 



C. Gluon Distribution in the DLA 

Similarly to the quark case, one can consider the vir- 
tual gg pair interacting with the nuclcons in the nuclear 
target. Using the well-known relationship between the 
gluon density and the deeply-inelastic cross section , 
the gluon density can be written as 



xgN{x, Q^) 



TT 



^ f dz f d^r\tP{z,r)\''<jggN{r) , (3.17) 



Jo 



where the wave function squared for a gluon pair is 
1 



|^(z,r)p = 



[{a^K2{ar) - aKi{ar)/rf 



z{l-z) 

WK,{arf/r^] , 
zr^ 

and the cross section in the DLA is (23j 
37r2 



(3.18) 



^ r^asiAy)x'g^^^ix'Ay). (3.19) 



Note (TggN — 9(TqqAr/4. Substituting (3.18) and (3.1£) 
into (3.17), one obtains 



9N{x,Q)=gN{x,Q,) + -J^ —4 

xa,(g'2)(^6^)5^i^^(x',g'2), (3.20) 

a form in which the small- a; limit splitting function in 
the DLA can be easily identified: Pgg{z) ~ 6/2: where 



z = x/x'. Comparing (3.2C) with ( [3.10| ), one finds that 
due to different splitting functions, the gluon density in- 
creases as increases 12 times faster than the quark 
density, leading to a muc h str onger scaling violation. 
Similar to the quark case, ( 3.20 ) is not to be interpreted 
as the integral equation since the right hand side does 
not contain f^^r but g^i"^ ■ Namely, g^^ is the gluon 
density seen by a gg pair with a frozen transverse sepa- 
ration while gf4 is the gluon density seen by the probe, 
e.g., the virtual Higgs boson. However, since the inte- 
gration kernel in ( 3.20 ) coincides with that of the DLA 
evolution equation, gjq should be very close to 

gDLA_ We 

confirm that numerically they differ only by few percent. 
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IV. GLAUBER MULTIPLE-SCATTERING 
MODEL 

The Glauber multiple-scattering theory treats a 
nuclear collision as a succession of collisions of the probe 
with individual nucleons within nucleus. A partonic sys- 
tem (/i), being a, qq ot a gg fluctuation, can scatter co- 
herently from several or all nucleons during its passage 
through the target nucleus. The interference between the 
multiple scattering amplitudes causes a reduction of the 
hA cross section compared to the naive scaling result of 
A times the respective hN cross section, the origin of 
the nuclear shadowing. In the high-energy eikonal limit, 
the phase shift of the elastic hA scattering amplitude is 
just the sum of phase shifts of individual hN scatter- 
ings. The total hA cross section is determined by the 
imaginary part of the forward elastic scattering ampli- 
tude as governed by the unitarity relation or the optical 
theorem. The unitarity of the total cross section is thus 
built-in within the Glauber theory. 

In the high-energy eikonal limit where the projectile's 
momentum k is much larger than the inverse of the inter- 
action radius, one describes the scattering amplitude by 
the Fourier transformation of the so-called profile func- 
tion r(b) in the impact parameter space 



(4.1) 



where k' is the outgoing wave momentum. In the high en- 
ergy scattering process, the collision is almost completely 
absorptive corresponding to an approximately real r(b) 
and pure imaginary amplitude. The total cross section is 
given by the optical theorem 



atot = Ylm./(0) = 2 Jd^hT{h). 



(4.2) 



To extend these results to nuclei, one assumes that each 
collision with the nucleon modifies the piece of the inci- 
dent wave passing through it by a factor 1 — Fjv (b — s) 
where s is the transverse position of the nucleon. The 
overall modification of the wave is thus given by 



n[i-r^(b-s,)] 



(4.3) 



Clearly, only those nucleons in the vicinity of impact pa- 
rameter b can influence Ta- The nucleus profile function 



is just the factor in ( |4.3| ) multiplied by the probability of 
finding nucleons in the positions ri, ...,ryi 



|^(ri,...,r^)|2=n^(^ 



(4.4) 



where p{Yi)(pYi is the probability of finding a nucleon at 
Tj. The nucleon number density is then pA(r) = P{^i)- 
For large A, the nucleus profile is calculated 



A 

F^(b) = 1 - n / d''v,p{v,)[l - Fjv(b - s,)] 
i=i •' 

A 

2 — 1 

J d^sTA{s)TN{h-s) 



(4.5) 



1 — exp 



where we have used A^ 1. We shall use a Gaussian pa- 
rameterization of the nucleus thickness function defined 
by 



T^(b) = / dzpAih,z) = 



7:R\ 



(4.6) 



where the mean nucleus radius R\ = QAR^g and 
Rws = 1.3^1/3 fm. 

Since the nucleon has a smaller transverse size than the 
typical nuclear dimension, one approximates the nucleon 
profile function in (4^) by Fjv(b — s) ~ ahN^^i^ — s)/2 
and obtains the total hA cross section 



a^A = / d^h2 



1 



,-ahNTA{b)/2 



(4.7) 



For a Gaussian parametrization of rA(b), the integration 
on the impact parameter can be done exactly leading to 



ah A = '2.7: R J 



1 2^ 



+ (_1)"-1^ 



2ttR\ + In Kh + Ei{Kh)] 



(4.8) 



where Kh = AahN /i2nR\), 7 = 0.5772 is Euler's con- 
stant and En{x) is the exponential integral function of 
rank n. The dimensionless quantity serves as an im- 
pact parameter averaged effective number of scatterings. 
For small value of k/i, ah a — > 2Tr R^Kh = AahN, the to- 
tal hA cross section is proportional to A or the nuclear 
volume. In the limit Kh 00, the destructive interfer- 
ence between multiple scattering amplitudes reduces the 
cross section to ah a 27ri?^(7 -f In k/i). Namely, the ef- 
fective number of scatterings is large and the total cross 
section approaches the geometric limit 27tR\, a surface 
term which varies roughly as A"^/^. In the kinematical 
region that we are interested in none of these limits are 
actually approached. The shadowing effect is determined 
by the strength of Kh , which is different for ah = qq than 
for a h = gg case. 

So far we have sketched the necessary assumptions 
for applying the Glauber theory for the elastic or the 
quasielastic scattering, that is, the longitudinal momen- 
tum remains the same during the passage of h through 
the nucleus. This condition can be relaxed to the for- 
ward production processes where there is a longitudinal 
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momentum transfer appropriate to the qq or the gg pair 
production. In this case, the scattering amphtude in (4.1) 
is modified to include the integration of the longitudinal 
position z due to the fact that an incident wave e'^^^^ 
for the photon 7 is changed to a final wave — r(b)e*'^^^ 
for particle h. In general, kj 7^ fc^ , and there is a phase 
difference between the 7 and h waves behind the target, 
i.e., e*^'^^^ where A/c^ = kJ — k'^. This leads to a coher- 
ence length cutoff which is governed by Ic — l/Akz- If 
Ic > Ra, the hadronic fluctuation h interacts coherently 
with all nucleons within the nucleus and the Glauber the- 
ory is valid. On the other hand, if Ic < Ra, h interacts 
coherently only with some of the nucleons and the inter- 
ference effect is reduced. 

To estimate the coherence length Ic, one calculates the 
change of the longitudinal momenta in the small-x region 



Afc, 



2v 



[x + x')m 



(4.9) 



where M is the invariant mass of the qq or the gg pair as 
defined in (3.7). Since the cross section is dominated by 
x' near x, the coherence length Ic — l/{2xm) exceeds the 
radius of heavy nuclei Ra ~ 4 fm as long asa;<2xl0~^ 
and the Glauber multiple scattering model is adequate. 
For a; > 2 X 10^^, the Glauber approximation breaks 
down and one should take into account the finite coher- 
ence length. Although we arc mainly concerned with the 
small-x behavior of the nuclear shadowing, we would like 
to get some handle on the coherence length cutoff in the 
moderately large- a; region. We propose to take a proba- 
bility approach [|l3|,|6) and introduce a coherence length 
form factor T^l^k^) = exp[-i?^(AA:^)2/4] ^ while pre- 
serving the nice structure of the Glauber formula. The 
probability for a double scattering is reduced by r, and 
that for a n-scattering is reduced by r"~^ so that one 



may make the following replacement in (4.7) 



1 



-o-hjvT4(b)/2 



1 



-o-h„TA(b)r/2 



(4.10) 



As T 1, the or igina l Glauber formula is preserved, 
while as r ^ 0, ( 4.10 ) gives ahA = AcrhN, i-e., there 
is no shadowing. In the kinematical region x < 10~^, 
the coherence length cutoff does not play a role. Even 
for X > 10~^, because the cross section is already very 
small, the effect is not very large. 



V. NUCLEAR QUARK DISTRIBUTION 



We now examine the effect of the gluon density in 



(3.12) on the shadowing of the nuclear structure function. 
Replacing ahN in (3.9) and (3.17) by ahA as calculated 
by the Glauber multiple scattering model, one obtains 
for (sea) quarks 



37?2 P fir' (^'^ 

xU{x,Q')^xfA{x,Ql) + ^x % dQ" 



X - [7 -I- In(Kgr) + El [KqT)] 



(5.1) 



where the effective number of scatterings k is given by 
2A-n 



Kg (X, Q ) 



3i?2Q2 



a,{Q')xg'ri'^^{x,Q') , 



(5.2) 



Expanding [• • •] in ( |5.l| ) and keeping only the first term 
which is proportional to Kg, one recovers the results for 
the nuclcon in the same DLA approximation. 

The nuclear structure function F^{x,Q'^) is defined 
t hro ugh the total ^*-A cross section, similar to in 
( |3.8|) . To improve the large- cc behavior, we add to 
the valence quark contributions and assume that there is 
no shadowing for the valence quarks 



F^{x, Q2) = J2 el[xfA{x, Q') + Axv,{x, 



(5.3) 



In principle, the initial condition xfA{x, Qq) in ( ^.l| ) may 
be calculated in the vector-meson-dominance (VMD) 
model. Fortunately, there have been some experimen- 
tal measurements on the nuclear quark shadowing in the 
low-Q^ region so that the initial condition can be ex- 
pressed in terms of the production of the initial nucleon 
structure function and the shadowing ratio Rq{x) 



Fi{x,Ql) ^ ARl{x)F^ {x^QD 



(5.4) 



where F^{x,Ql) is given in ( |3l^ ). Recently E665 Col- 
laboration has measured the nuclear shadowing effect in 
muon scatterings off ''"Ca and ^°*Pb at low a; B. In 
general, the small x data are measured at low Q , e.g., 
for a; = 6 X 10"'' the average (Q^) = 0.26 GeV^ while 
a: = 6 X 10^^ corresponds to (Q^) = 1.35 GeV^. We 
parametrize the initial shadowing ratio Rq{x) at Qg = 0.4 
GeV^ and use (5.2) to calculate the nuclear structure 
function at the measured (Q^) values at different x val- 
ues. The nuclear shadowing ratio is defined by 



i?,(a:,Q2) 



AFi'ix,Q^) 
i?0(a;)Ff(x,Q2) 



AF,^{x-Q^,Ql) 



F2 {X, Qq) 



AF^{x; 



\Ql) 



(5.5) 



where AF2 = F2 ~ F2{Qq) can be explicitly calculated 
by performing the double integration over x' and Q'^ in 
T1). The results for "^^Ca and ^°^Pb are shown in Fig. 

The lower panels of Fig. ^ show our parametrizations 
of initial shadowing ratios at Qg evolved to the measured 
(Q^) compared to the experimental data from E665. For 
(Q^) < Qq, we take it to be at the minimal scale Qg. 
As is evident from the comparison, our parameteriza- 
tions faithfully represent the experimental data at small 
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X. The shadowing ratios Rq{x, Q^) at high values of 
are shown in the upper panels of Fig. ^. 
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M E665 '95 

- Evolved to <qV, 
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FIG. 4. The a;-dependence of nuclear quark shadowing ra- 
tio Rq = F.f/iAFi') at various for ^°^Ph and '"'Ca. The 
lower panels represent the parameterizations of initial shad- 
owing ratios evolved to the experimental (Q^)exp values. The 
non-saturation of Rq{x) as x decreases at large is the fea- 
ture of the perturbative shadowing mechanism. 

We note that at large the nuclear shadowing ef- 
fect is reduced but does not diminish, i.e. it is not a 
higher twist effect. This can be understood as the inter- 
play between the perturbative and the non-perturbative 
shadowing mechanisms. The nuclear quark density is 
initially shadowed by i?° at Qq compared to the free nu- 
cleon case due to some non-perturbative mechanism such 
as coherent scatterings of the p-meson as discussed in the 
VMD model As increases from Qq, the nucleon 
structure function increases by AF^ {x;Q'^,Qq) as gov- 
erned by the GLAP evolution equation in the DLA. If 
the nuclear increases by the same amount (times A) , 
i.e., we neglect the multiple scattering effect and evolve it 
according to th e sam e GLAP equation, the initial shad- 
owing effect in (5^) decreases quickly due to the strong 



scaling violation at low {AF^ {x; Q^, Qq) is not small 
compared to F^{x,Qq)). This is illustrated in Fig. |^ 
where the shadowing ratio is calculated by taking the 
first term in the expansion [• • •] in (^^) and is plotted as 
a function of a.t x = 10"''. The shadowing effect dis- 
appears quickly when both nucleon and nuclear structure 
functions evolve according to the GLAP equation. 

However, because xg^^^ increases rapidly as in- 
creases, the growth of the gluon density in a nucleon 
compensates for the suppression due to the smallness of 



the size of the qq pair. As a result, Kg is not small in the 
semihard region. In fact. Kg increases as increases 
from Qq to 2 GeV^ and becomes greater than 1. The 
perturbative rescattering effect becomes important and 
needs to be taken into account in the Glauber formal- 
ism. The effect of multiple scatterings slows down the 
Q^ evolution of F^ as compared to F^ , i.e.. 



AF^{x;Q',Q'q) < AF^{x-Q\Ql) , 



(5.6) 



which provides a perturbative shadowing mechanism be- 
yond Qq. This is illustrated in Fig. |^ where the Q^ evo- 
lution is indeed slower than the GLAP evolution case. 
As Q^ becomes even larger (Q^ > 5 GeV^), Kg decreases 
fast and the probability for rescattering becomes small. 
In this case, both F^ and F^ evolve slowly and so does 
the shadowing ratio. Also plotted in Fig. |^ is the double 
scattering case where we take the first two terms in the 
expansion 



[7 -I- luK, -I- El{Kq 



4N 



(5.7) 



As it is evident from Fig. ||, the double scattering over- 
estimates the amount of shadowing. Nevertheless, it ac- 
counts for most of the multiple scattering effect. 
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FIG. 5. The Q^-dependence of nuclear quark shadowing 



ratio Rq = Fi/{AF^ 



at s = 10-" for ^''^^Pb 



The interplay between soft physics and the partonic 
picture of the nuclear shadowing is also evident in the 
x-dependence of Rq in Fig. ^. The apparent flatness of 
the shadowing ratio at low Q^ in the small- a; region is 
altered by the perturbative evolution. This is due to the 
singular behavior of xg^^ as a; ^ at large Q^ leading 
to the strong a;-dependence of the effective number of 
scatterings, Kq{x). Thus, the saturation of the Rq as x 
decreases is a feature only at low (< 2 GeV^) where 
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the gluon density is not so singular, and we do not expect 
saturation to happen for higher values of . 



VI. PREDICTING GLUON SHADOWING 

In the previous section, we have emphasized the im- 
portance of the interplay between the non-perturbative 
and perturbative shadowing mechanisms in order to pre- 
dict the nuclear quark shadowing at large . This would 
suggest that the prediction for the nuclear gluon shadow- 
ing is uncertain because there is very little experimental 
information on the gluon shadowing at low . However, 
as we shall demonstrate below, this is not the case and 
the gluon shadowing can be unambiguously predicted for 
large in the context of perturbative QCD. 

The nuclear gluon distribution can be calculated in the 
Glauber model similarly to the quark case p^ ] 



x9a{x,QI) -f 
1 



Jx JQl 



X - [7 -I- In(KgT) -t- El{KgT)] 

T 



(6.1) 



The main difference between the quark and gluon cases 
is different wave functions of the virtual probe and the 
effective number of scatterings, Kg — 9Kq/4 due to dif- 
ferent color representations that the quark and the gluon 
belong to. As a result, the gluon density increases as 
increases 12 times faster than the sea-quark density at 
small X. The two important effects which make gluon 
shadowing quite different from quark shadowing are a 
stronger scaling violation in the semihard scale region 
and a larger perturbative shadowing effect. This can be 
seen by considering the shadowing ratio 



Rgix,Q^) 



xgAix,Q'^) 
AxgNix, Q2) 
x9Nix,Ql)R°{x) 



(6.2) 



xgN{x, Ql) -f AxgN{x; Q^,Ql) 



where Rg{x) is the initial shadowing 
Axg{x;Q'^ ,Ql) is the change of the g' 
as the scale changes from Qq to Q^. 
ing violation due to a larger Kg at 
AxgNix;Q'^,Ql) > xgN{x,Q'^) as 
than 1 2 GeV^. The dependence of 
initial condition R^{x) diminishes and 
shadowing mechanism takes over, i.e. 



ratio at Qq and 
luon distribution 
The strong seal- 
small X causes 
becomes greater 
Rg{x, Q^) on the 
the perturbative 



Rg{x,Q^) 



AxgA{x;Q'^,Ql) 
AxgN{x]Q'^,Ql) 



(6.3) 



for Q2 > 2 GeV2. 

The above statement can be made more rigorous by 



taking a partial derivative of (6.2) with respect to InQ'^ 
at fixed x 



^^^^^ = 7(x,g^)[<T(a;,Q^) 



91nQ2 



-Rg{x,Q^ , 



(6.4) 



where the anomalous dimension 7(0;, Q^) and the so- 
called perturbative shadowing ratio i?g are defined 



7(2;, 



91n xgN{x, Q^) 
ainQ2 

dxgA{x,Q'')ld\nQ'' 
AdxgN{x,Q^)/d\nQ'^ 



(6.5) 



Note that Rg is the ratio of the derivative of xgA to 
that of xgN , which is independent of the initial shadowing 
ratio and can be calculated perturbatively. In the small- a; 
region 



R. 



PT 



7 + lUKg + El{Kg) 



= 1 _ ^ + ^ 

4 18 



(6.6) 



The interesting feature of (6.4) is the negative sign in 



front of Rg on the right hand side such that it is self- 
stabilizing. If i?g < at g2 and dRg/dlnQ'^ > 0, the 
shadowing ratio increases with Q^; while for Rg < R^'^ , 
the shadowing ratio decreases with Q^. As a result, the 
Q^-evolution of Rg is such that it is always in the process 
of approaching R^^ ■ Because R^'^ also changes with Q^, 
the Q^-evolution of Rg never stops. The anomalous di- 
mension 7(x, Q^) determines how fast Rg is approaching 
R^. 
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FIG. 6. The Q -dependence of the anomalous dimension 



7(a;,Q2- 
X = 10" 



and the perturbative shadowing ratio R^"^ {x,Q^) 



at 



for Pb. Also plotted is the shadowing ratio Rg 



with an initial condition iig = 0.5. 



Introduce an important scale Qgjj = 3 GeV^ ("SH" 
stands for "semihard") where j{x,Qg^) ~0(1) for small 
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values of X (10 



<x<10 ^). In the semihard region 
7(x, Q^) is large and Rg approaches 



i?g ^ (x, Qgjj) fast regardless of what the initial shadow- 
ing is at Qn. It is precisely in this region that k 



IS 



large and (^^) gives a large shadowing correction. As 
> QsH' anomalous dimension 7 gets small and 
in the meantime R^"^ approaches 1 rather quickly, the 
shadowing ratio Rg evolves slowly. Thus, the behavior of 
the Q^-evolution of the gluon shadowing can be under- 
stood by the interplay between the anomalous dimension 
of the nucleon gluon structure function and the pertur- 
bative shadowing mechanism. 

In Fig. ^ we show the Q^-behavior of 7(0;, Q^) and 
i?P^(x,Q2) at X = lO"'^ for ^ospb. The evolution of the 
shadowing ratio Rg with an initial condition -R°(a;) = 
0.5 i s also plotted to illustrate the evolution governed by 

To study the sensitivity of the shadowing ratio on the 
initial condition, one can solve ( |6.4|) explicitly by making 
a transformation Rg = R'gexp[— J ^(x,t')dt'] where t — 
InQ^. Then R'g can be directly integrated out which 
yields 



R, 



,{x,t) ^ Rg{x,tn)e ■''0 +e ■''o (6.7) 



X / 7 

'to 



)dt' 



In the absence of perturbative shadowing, i.e. R^"^ — 1, 



the integral in { p.7\) can be done exactly leading to a pure 
GLAP evolution 



R 



GLAP 



(x,t) = 1 - [1 - Rg{x,to)]e 



- r* i(x,t')dt' 

J tn 



In the DLA, the integration of 7(0;, t) yields 



^^^ix,t) My) 



(6.9) 



where y is defined in ( ^.13 ). The factor l/Io{y) sup- 
presses the initial shadowing Rg{x,to) quickly for large 
t > tsH = InQgjj- The integration in (6.9) picks up 
most important contributions in the region to < t < igH- 
1 quickly as t > igu- In the pres- 



Therefore, i^^^AP 



bigger as x gets larger due to a weaker scaling violation, 
which is evident in Fig. |^. 
1.0- 




Q' (GeV') Q' (GeV') 

FIG. 7. The maximal uncertainty in predicting the gluon 
shadowing due to initial conditions. The shaded regions are 
bounded by two extreme initial conditions: R'^{x) — 1 and 
i?«(x) =0. 



Neglecting the first term in (3/7) and approximating 



the perturbative shadowing by R^^ = R^'^{Qsh) 



Ql<Q^ < QIh and i?^^ = 1 for q2 > 



one obtains 



from (|6.7|) the shadowing ratio for Q > Qgjj 



Rg{x,t) = 1- [1-i?^ 



PT 



(X,tsii)\e -"sn . (6. 



10) 



Comp aring ( 6.IOI ) with (6.7), one notices that solution 



(|6.10| ) is just solution (6.8) with a new initial condition 



ence of perturbative shadowing, i.e. Rg < 1, Rg{x,t) 

at large t is mainly determined by the second term in though R^J{x, Q^) approaches 1 quickly as 



i?o ^ R^g^{x,Qly) and Ql replaced by The dif- 

ference is that the evolution factor exp[— /] no longer 
falls rapidly and Rg evolves rather slowly. Thus, the Q^- 
behavior of gluon shadowing at large can be roughly 
described by an initial condition at Qg^ followed by a 
GLAP evolution. The important point is that the initial 
condition at = OIh is pre-determined by the per- 
turbative shadowing. Our main conclusion is that even 



( |6.7[ ) and it depends only on 7 and -Rg"^, especially in 
the small- a; region. 

Fig. 1^ demonstrates the stabilization of the gluon shad- 
owing for any initial conditions as increases at various 
values of x for ^°^Pb. The shaded regions indicate the 
maximal uncertainty in predicting the gluon shadowing 
as they are bounded by two extreme initial conditions: 
Rg{x) — 1 and Rg{x) — for all values of x. One can 
almost uniquely determine the amount of the gluon shad- 
owing for > 3 GeV^. In general, the uncertainty gets 



increases 

beyond Qgjj (a higher twist effect), the shadowing ratio 
at large is still governed by the perturbative shadow- 
ing at the semihard scale, namely i?^"'"(x, Qgjj). 

The x-dependence of the gluon shadowing can also be 
predicted as long as > Qgjj where the influence of 
the initial condition is minimal. Fig. ^ shows the x- 
dependence of the gluon shadowing ratio at various 
values. The shaded region is bounded by the distribu- 
tions calculated using two drastically different initial con- 
ditions Rg{x) — and Rg{x) — 1. It is clear that the 
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shape of the distribution is quite robust in the small- 
X region regardless of what initial conditions one may 
choose. Due to the perturbative nature of the shadow- 
ing, these distributions do not exhibit a saturation as x 
decreases. This is related to the singular behavior of the 
nucleon gluon distribution at large and small x. As 
long as the nucleon gluon density does not saturates at 
small X, the shadowing ratio does not saturate either. 
This is in accordance with the observation by Eskola, 
Qiu and Wang in the context of the gluon recombination 
model [||. The rapid Q^-evolution of the gluon shad- 
owing ratio in the low region found by Eskola p^ ] 
can be understood in our approach as the quick onset of 
the perturbative shadowing. In Eskola's approach some 
shadowing persists even at = 100 GeV^ which can 
also be understood in our picture as the region of fully 
developed GLAP evolution. We note that multiple scat- 
terings inside the nucleon are also possible in the region 
of very small x, where the eikonalization of the nucleon 
cross section is needed. However, in this region, the DLA 
breaks down and a BFKL resummation is necessary. 



^.(y 



0.9- 



X 

< 



: 0.8- 



•0.7- 



0.6- 



0.&^ 







Q' (GeV') 

















10 



10^ 



10^' 



10 



FIG. 8. The a;-dependence of the gluon shadowing ratio at 
various values. The shaded region indicates the uncer- 
tainty due to the initial conditions. 



on central events where the impact parameter is small 
(e.g. 6 < 2 fm). It is thus necessary to study the impact 
parameter dependence of the nuclear shadowing effect. 

Intuitively, one expects the nuclear shadowing effect to 
depend on the local nuclear thickness at a given impact 
parameter. In the Glauber approach, the effective num- 
ber of scatterings is the product of the cross section and 
the nuclear thickness function. Since the shadowing is 
a non-linear effect in the effective number of scatterings, 
the impact parameter dependent shadowing ratio can- 
not be factorized into a product of an average shadowing 
ratio and the nuclear thickness function WA- 



1.0 
0.9 



DC 



-q 0.8 - Average 
q . /x=io- 

4 0.7 \ / Q'=10GeV' 

R. 

0.6 

Q 5 I \ \ \ \ \ \ \ \ L_il \ \ \ L_ 

■ 1 2 3 4 5 6 7 
Impact Parameter b (fm) 




Average 
b=6fm - 

- - b=3fm 

- - - b=0 fm ' 



°Pb 
Q'=10GeV' 



FIG. 9. The impact parameter dependent gluon shadowing 
ratio calculated assuming Ul(x,Ql-h) = 1 for any h. The 
upper panel is the shadowing ratio as function of 6 at a; = 10~* 
and = 10 GeV'^. Tie lower panel shows the a;-dependence 
of the shadowing ratio at — 10 GeV^ and various values 
of 6. 



VII. GLUON SHADOWING AT FIXED IMPACT 
PARAMETER 

So far we have confined ourselves to the impact param- 
eter averaged nuclear structure function. However, in the 
binary approximation of hard scattering processes in pA 
and AA collisions, one has to calculate the parton cross 
section at the nucleon-nucleon level at each impact pa- 
rameter. In addition, to enhance the possible signal of a 
quark-gluon plasma in heavy-ion collisions, one triggers 



The impact parameter dependent nuclear structure 
hmction xgA(x^Q'^\\i) is defined as 



(x,Q2) = y d^hx~gA{x,Q'\^) 



x2 1 



,-CrggJvTA(b)/2 



(7.1) 



Comparing (7T) with (4/7), one finds 

4 dx' ['^^ 
xgA{x,Q^;h) = xgA{x,Ql:h) + — — dQ'^ 

Jx Jq2 

(7.2) 
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where xgA{x,Qo',i>) is the initial condition at fixed im- 
pact parameter b. As shown in the previous section, the 
prediction for gluon shadowing ratio at > 3GeV^ does 
not depend on the initial condition. We take the initial 
shadowing ratio to be 1 for all b's, i.e. no initial shadow- 
ing. The impact parameter dependent shadowing ratio 
is thus defined 



TA{h)xgNix,Q^) 



(7.3) 



In Fig. ^ we show the gluon shadowing at various im- 
pact parameters. The upper panel is the shadowing ratio 
as function of 6 at a; = lO"* and ^ iq GeV'\ The 
parton density is more suppressed at a small impact pa- 
rameter. The x-dependence of the impact parameter de- 
pendent shadowing ratio is also shown in the lower panel 
of Fig. I 



VIII. CONCLUSIONS 

We have studied the nuclear shadowing mechanism 
in the context of perturbative QCD and the Glauber 
multiple scattering model. The nuclear shadowing phe- 
nomenon is a consequence of the coherent parton multi- 
ple scatterings. While the quark density shadowing arises 
from an interplay between the "soft" and the semihard 
QCD processes, the gluon shadowing is largely driven by 
a perturbative shadowing mechanism due to the strong 
scaling violation in the small- a; region. The gluon shad- 
owing is thus a robust phenomenon at large and can 
be unambiguously predicted by perturbative QCD. 

The scale dependence of the nucleon structure func- 
tion has been used to put a stringent constraint on the 
gluon density inside the nucleon in the double logarithm 
approximation. It is shown that the strong scaling vio- 
lation of the nucleon structure function in the semihard 
momentum transfer region 0.8 < < 3 GeV^ at small 
X can be reliably described by perturbative QCD and is a 
central key to the understanding of the scale dependence 
of the nuclear shadowing effect. 

We have also studied the impact parameter depen- 
dence of gluon shadowing and shown that it is a non- 
linear effect in the nuclear thickness function. It is im- 
portant to correctly incorporate the impact parameter 
dependence of the nuclear structure function when one 
calculates the QCD processes such as minijet production 
and J/i/j production in the central nucleus collisions. 
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APPENDIX: PARAMETRIZATION OF GLUON 
SHADOWING AND ITS IMPACT 
PARAMETER DEPENDENCE 

Motivated by our studies of the gluon shadowing at > 
QsH' parametrize the gluon shadowing ratio in the 
following way: 



Rg{x,Q^) = l-[l-R^g^{x,Ql^)] 
X [l-aolnQVQ2jj] 

where Qgjj = 3 GeV^, gq = 0.3 and 



where A is the atomic number of the nucleus and 
ai = 0.0249 ; a = 0.18 ; A = 0.162 ; /3 = 18 



(1) 
(2) 
(3) 



The impact parameter dependent ratio Rg{x, Q^; b) is 
parametrized via a 6-averaged ratio Rg{x, Q^) 



Rgix, Q';h) = 1 - co(l + ciA-" b^) 



(4) 



X exp[-A-"'&Vc2](l - Rg{x, Q^)) , 

where h is the impact parameter in fm and 

Co = 1.7 ; ci = 0.043 ; a' = 0.11 ; ca = 11.45 . (5) 

Note that these parametrizations are valid for > Qgjj 
and X < 10^^. 
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